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A semi-phenoinenological approach to describe the evolution of Fermi 
surface (FS) and electronic structure with doping is presented which is based 
on the spin-ferniion model. The doping is simulated by a frustration term in 
the spin Hamiltonian and the complex internal structure of the spin-polaron 
quasiparticle is taken into account by a superposition of spin-polaron states 
with different radii. By calculating the spectrum, the spectral weights and 
the FS we find a rather drastic change of the electronic structure with 
doping which explains many photoemission data, i.e. the isotropic band 
bottom and the remnant FS of undoped cuprates, the large FS and the 
extended saddle point of optimally doped compounds and the pseudogap 
'^ . in the underdoped samples. 

§ : PACS number(s): 75.50.Ee, 74.20.Mn, 71.38.+i, 75.30.Mb 

^ : I. INTRODUCTION 

^ ■ The evolution of Fermi surface (FS) and electronic spectrum in the normal state of high 

Qv^ , temperature superconductors are intensively studied at present. Recent angle- resolved 

O I photoemission spectroscopy (ARPES) studies indicate different dispersion relations in the 

insulating and optimally doped cuprates. The undoped compounds [jl],|| show an isotropic 

band bottom close to A^ = (7r/2, 7r/2) in momentum space, a large energy difference 

between A^ and X = {it, 0), and a so-called remnant Fermi surface (that surface where 

the single particle spectral weight shows a sudden drop of intensity). In the optimally 

^ [ doped compounds a flat band region, a large FS centered at M = (vr, vr) and a so-called 

5 ' shadow Fermi surface due to antiferromagnetic (AFM) spin correlations |@ are found 

O ■ |]5Hri|. The fiat band region has the form of an extended saddle point which stretches 
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in the direction (7r/2, 0)-(7r, 0). The shadow FS resembles the main FS but is shifted by 
'k><( I the AFM vector Q = (7r,7r). Furthermore, the underdoped compounds demonstrate the 
^ ' existence of a high energy pseudogap near the point X with an energy of about 0.1 -t- 0.2 
eV PJT^-|l^ • Assuming a simple rigid band filling, the isotropic minimum of the undoped 
compounds would suggest small hole pockets around A^. However, there are no clear 
experimental indications for them [|H| and the Fermi surface seems to develop in an arc- 



like fashion |T^. It is a real theoretical problem to reconcile that seeming contradiction 
which can be attacked only by considering the spectral density evolution as a whole. 



There are numerous studies of the hole spectrum in two dimensional (2D) doped 
antiferromagnets using a wealth of theoretical models, like the generalized t-J, the effective 
three-band model, the Kondo-lattice or the Hubbard model. Most of those studies are 
restricted to a very low hole concentration using the exact diagonalization of small clusters 
|lT6| , p!7| , the quantum Monte Carlo technique |1^, the self-consistent Born approximation 



(SCBA) [p!9| -|2T[| or a "string" ansatz for the hole wave function |2^. One can start either 
from the two-sublattice Neel-type state [|l^,^,^,^ or from the spin rotationally invariant 



spin liquid state [^Ijil] and obtain qualitatively the same result: the hole motion takes 
place mainly on one sublattice, i.e. the dispersion relation is dominated by an effective 
hopping to the next nearest neighbors and the dispersion has its minimum at A^ for small 
doping. There are much less attempts to explain the Fermi surface evolution, like for 
instance p4| which used a diagrammatic method combined with a semi-phenomenological 
spin susceptibility. The exact diagonalization studies (on clusters containing up to 16 -=-20 
sites) at finite hole doping [^5H 27[| have difficulties to reconstruct the FS and to investigate 



a sufficiently smooth variation of doping due to the restricted set of k-points. 

The approach which shall be presented here is semi-phenomenological as well. On one 
hand, we start with the spin-fermion model (the reduced form of the three-band Emery 
model) that has a clear microscopic basis [^,|^. We also include direct oxygen-oxygen 
hopping terms which are rather important. On the other hand, we simulate the doping by 
a frustration term J2 in addition to the nearest neighbor exchange Ji of the spin part in 
the Hamiltonian. The similarity between doping and frustration had been ffist proposed 
in [^] based on a similar decrease of the magnetic correlation length. It is interesting to 
note that cluster calculations indicated a considerable value of the frustration parameter 
J2I Ji ~ 0.1 even for the undoped La2Cu04 ||3l|. Of course there is no full equivalence 
between doping and frustration. For example, the doped t-J model and the frustrated 
J1-J2 model give different results for the dynamical spin-spin structure factor and for 
the spectrum of Raman scattering [^. But numerical calculations on finite lattices 
indicate the equivalence of the mentioned models if we are interested in the static spin- 



spin correlation functions [Q . And those quantities are especially relevant in our present 
Green's function method where the spectrum £:(k) is determined by the static spin-spin 
correlation functions of the spin subsystem. 



Previous studies of the frustrated generalized t-J [|J] and the frustrated effective 
three-band model [^ demonstrated already the strong influence of frustration on the 
spin-polaron spectrum. However, these investigations were performed only within the 
local polaron approach (LPA) which cannot describe the influence of long range spin order 
on the excitation spectrum. Furthermore, only the bare spectral weight was considered 
which could only give a hypothetical FS. Recently [^ the frustrated t — t' — t" — Ji — J2 



model was discussed in the framework of the SCBA where the holes are treated as spinless 
fermions and the spins as normal bosons I^Oj. The resulting SCBA in [^ corresponds 



to a two sublattice approach in which the spectrum is dictated by the symmetry of the 
reduced magnetic Brillouin zone (BZ) (e.g. the points F and M are equivalent). However, 
such a symmetry is not reproduced in the ARPES data. On the contrary, the spherically 
symmetric approach (in the spin space) is used in the present investigation to describe 
the spin subsystem. As a result, the spectrum is periodic relative to the full Brillouin 



zone (BZ) of the Cu02 plane. 

The distinctive feature of our present investigation consists in treating the spin-polaron 
as a complex quasiparticle, i.e. as a superposition of spin polaron states with different 
radii |^. Our basis operators enclose the local polaron (bare holes and the Zhang- 



Rice polaron) and also its coupling to spin waves with quasimomentum q close to the 
antiferromagnetic wave vector Q. By introducing the complex structure of the spin- 
polaron one may describe the important splitting of the lowest local polaron band. This 
allows, for example, to reproduce for heavily underdoped compounds a sudden drop in the 
intensity of ARPES peaks as k goes from A^ to M. And we will show that it reproduces 
also many of the other features of electronic structure evolution mentioned above. 



II. THE HAMILTONIAN AND THE METHOD 

The main features of the hole motion in the Cu02 plane are described by the model 

H = f + J + h, (1) 



28,29 



f = 4r^pU- + 5RUR 



R 



J - ^Yl SaSR+g + ^Yl SRSR+d (3) 

R,g ^ R,d 

where Ji, J2> are the antiferromagnetic interactions between the first-nearest neighbor 
(g = ±g^, ±gy) and the second-nearest neighbor (d = ±g^±gy) on a square lattice. The 
direct oxygen-oxygen hopping is included by 

h=-hJ2 [4+a. (cR+a, + CR_a, + CR+g^+a, + CR+g^-a,) + h.c] , (4) 

R 

with 

PR = o E CR+a, Sn = 5^<5-", {pr, p^,} = (5r,r' + jJ2 '^R,R'+g ' (5) 

a ^ g 

where a.:r:,y = lgx,y 

Here and below the summation over repeated indexes is understood everywhere; 
{...,...}, [...,...] stand for anticommutator and commutator respectively; gx,y are basis 
vectors of a copper square lattice ( |g| = 1), R+a are four vectors of O sites nearest to the 
Cu site R; the operator cj^+a creates a hole predominantly at the O site (the spin index 
is dropped in order to simplify the notations); d" are the Pauli matrices; the operator S 
represents the localized spin on the copper site. We do not introduce the explicit relative 



phases of p- and d -orbitals since they can be transformed out by redefining the operators 
with phase factors exp (^QR) . In order to compare our results with other authors and 
with experiment we restore these phases by changing in the final results k — i> k' = k — Q. 
The parameter r is the hopping amplitude of oxygen holes that takes into account the cou- 
pling of the hole motion with the copper spin subsystem. For convenience, the exchange 
interactions Ji and J2 can be expressed in terms of the frustration parameter p: 

J^ = {l-p)J, J2=pJ, 0<p<l, J >0 . (6) 

The frustration parameter p has to be connected with the number x of holes per copper 
atom. We are guided by an estimate based on the single-band Hubbard model with 
realistic values of U/t ~ 5 for x = 0.1 where one finds a value oi p ~ 0.1. Let us note 
that in the case of La2_a;Sr2,Cu04 the spin system of the CUO2 plane looses its long-range 
order for x > 0.02. 

The adopted effective spin-fermion Hamiltonian can be obtained by canonical trans- 



formation from the Emery model |3^ in the regime Ud '^ t ':^ e, e = Ep — e^ > 0, 
where Ud is the Coulomb repulsion on copper sites. The values of the energy parameter 
set e, t, Ud niay be extracted from band structure or cluster calculations. Two typi- 



cal sets are: e = Ep — Ed = 3.6 eV, t = 1.3 eV, Ud = 10.5 eV, |]39[, or e = 2.0 eV 



t = 1.0 eV Ud = 8.0 eV iO|,|Tl. Taking into account that r ~ t'^/E, J ~ 4r(t/£ 



,2 



we 



adopt the parameter values r ~ 0.4 eV, J = 0.4r and h = 0.4r. Below we take r 
as the unit of energy. We consider a spin-liquid state with spin rotational symmetry 
for the copper subsystem, e.g. the spin-spin correlation functions satisfy the relation 

r< — / ca oa \ _ Q /c^fe'^) C^(?/'^)\ ^ Ca ^ n 

^r — \^R^R+r/ " ^ \^R ^R+r / 5 ^ ^R ^— ^■ 

To treat the problem of a spin-polaron in the projection method approach we introduce 
a finite set of basis operators connected with each cell R: A^ ^ (with spin a) , where i is 
the number of the operator, i < N. The set contains the creation of a hole (with spin 
a) and other operators which describe a hole connected with spin excitations at different 
distances. The choice of the set is dictated by the physical sense of the problem and will 
be explained in detail below. To calculate the spectral function of bare holes we take the 
sum of two contributions Cj^+a^ and c^^+a ^^at are related with two possible positions of 
an oxygen hole in the unit cell. 

As usual we introduce the retarded two-time Green's functions Gij{t,'k) defined in 
terms of the Fourier transforms Ak,! of the operators Ajif. 

G.,(t,k) = (A^MKji^)) = -^©W ({^k.W >^k,(0)}) , (7) 

viV R 
The equations of motion for the Fourier transforms of the Green's functions have the form 
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K,,{k) = {{A^^, , Al^}) , B^ 



In the projection technique we approximate the new operators -Bk,i by their projections 
on the space {Ak,i} of the basis operators: 

B^,^Y.Lu{^)^^,i, L(k) = D(k)ir-^ A,(k) = ({5k,,^,^,,}). (9) 

I 

After substitution of the approximate expressions for the operators i?k,i in (0) into the 
equations of motion (H), the equation system (j^) for Green's functions (^k,il^kj) be- 
comes closed and can be presented in the matrix form 

(loE - DK-^) G = K , (10) 

where E is the unit matrix. The quasiparticle spectrum e{h) is determined by the poles 
of the Green's function G and can be derived from the condition 

detlKeik) - D\ = . 



III. THE BASIS OPERATORS AND APPROXIMATIONS 

At first we choose a set of three basis operators which describe local spin-polaron 
states 

^R,l = ^{c-R+a, + CR-aJ ; ^R,2 = -^{cR+ay + CR_a,J ; ^R,3 = ^R^R . (11) 

The combinations of them constitute the Zhang-Rice polaron and bare hole states, in 
particular: 

A A \~^ 2kR 

(12) 

The spectrum of elementary excitations in the framework of the local polaron approach 
(LPA) was investigated previously [^. In particular it was pointed out that the frus- 
tration in spin subsystem and oxygen-oxygen hopping can explain the appearance of an 
extended saddle-point in the spectrum. The main shortcoming of LPA consists in the fact 
that the spectrum of elementary excitations depends only on the short range spin-spin 
correlation functions even at zero temperature. This means that LPA does not describe 
the influence of spin long range order on the excitation spectrum. 

Treating the periodic Kondo problem on a square lattice in the two sublattice spin 
structure, Schrieffer |^^ emphasized the crucial importance of taking into account coher- 



ence factors related to the existence of long range order. In this model the band electrons 
(holes) are coupled to localized spins and the simplest Hamiltonian has the form 



Hk = E^g4+gCR + /E4^RCR + ^^E^R+g^R (13) 

R,g R '^ R,g 

where the insite exchange with constant I is analogous to the f Hamihonian (0) in our 



model. In Ref. |^2| the mean-field approach was taken for the Neel state at T = and 



the spins were treated as classical vectors: 

5r = (5„,Soe*Q^, ^0 = const. (14) 

In this approximation the Kondo-interaction Hamiltonian (|13]) takes the form of a poten- 
tial energy with doubled period. As a result, this interaction leads to the hybridization 
of bare particle states with momenta k and k + Q and one should perform a standard 
u — V transformation to take into account this hybridization from the very beginning. In 
the adopted Neel state the amplitude Sq of the spin wave with q = Q (the Q-wave) 
has a macroscopic large value and has properties analogous to the amplitude of a Bose 
particle with zero momentum in the superfiuid Bose-gas. As a result, this amplitude can 
be treated as a c— number for many problems. Then the hybridization corresponds to 
coupling of the Q-wave to local electron states. But it does not represent new states and 
leads only to mixing of states with momenta k and k + Q. 

The distinctive feature of the present investigation consists in considering the one- 
particle motion on the background of the spherically symmetric spin subsystem state. In 
this background the average value < Sq >= and the above mentioned treatment fails. 
At T = and without frustration the only value which can be treated as a macroscopic one 
is < SqSq >. Then the coupling of a local state to Sq corresponds to a new delocalized 
spin polaron states - local polaron states (pUD dressed by antiferromagnetic spin wave Sq. 
For the model under discussion at T = such states were introduced in [^] and have the 
form 

QrA^y, 1 = 1,2,2, (15) 

Qj, = e-^Q^5Q ; 5q = S'^a^ ; Sq = N'^ E ^^'^''^Sr, . (16) 

Ri 



It was shown in |4^ that it is important to take into account the quantum nature of the 
spin Q-wave because the transitions between the local polaron states and the delocalized 
polaron states lead to a splitting of the lowest LPA bands and change their properties in an 
essential way. Only the complex spin-polaron approach, i.e. that approach containing the 
local spin-polaron operators plus those operators dressed by the Q-wave, can be compared 
with a Green's function calculation for the bare hole spectral function ^^(k, cu) taking 
into account the spin-polaron damping by means of the SCBA [Q. One can see that 



the SCBA quasiparticle peak and its intensity are well given by the lowest band of the 
complex spin-polaron approach, whereas the excited bands refiect the non-coherent part 
of A,,(k,a;). 

In the present paper we treat the model at nonzero temperatures T and frustration 
p. In this case, the spin subsystem looses its long-range order and the average value 

6 



< SqSq > is zero. So one must introduce new operators in order to preserve the complex 
spin-polaron results and to take into account a finite spin correlation length ^. It is natural 
to generalize the above choice of operators by introducing local polaron states dressed by 
spin waves with momenta q close to the antiferromagnetic vector Q (for such q the spin- 
spin structure factor is sharply peaked even at T, p ^ 0). These operators can be taken 
in the form: 

^R,4 = Qb. ^R,i; ^R,5 = Qr. ^r,2; ^R,6 = Qn ^R,3 (17) 

Qn = N-' E e'^^'Sn+p-, fl = {q, | ± tt - q,J < Ko} (18) 

Here fl is the square region around Q and around equivalent points (described by the 
four squares Q = kq x kq related to the corners of the first Brillouin zone (BZ), see Fig. 
1.1a). The choice of the parameter kq is dictated by the spin correlation length ^ which 
depends on the frustration. In order to clarify our future choice of the parameter Q{k,o) 
for different frustrations let us represent each of the operators in ( p!?] ) (taking A-r^^ as an 
example) in the form 



^R,4 = E «(R - ROe^'^^^-^') ^r'Ar,i (19) 

R' 

where 

The absolute value of a(R — R') depends only on the absolute value of the difference 
1 = R — R' and decreases with the increase of /. This dependence a(I) = |a(l)/a(0)| has 
the form 

a(I) = \sm{l^Ko)sm{lyKo)/{lxly^)\ (21) 

The value Iq , which satisfies the condition 

IqHq = 71, (22) 

can be qualitatively treated as a coupling radius of the local polaron which defines the 
range of coupling to the long distance spin correlations, since a(l) ~ 1 for / < Iq and 
(5(1) -C 1 for / :» /o- On the other hand, the coupling radius of a spin polaron must be 
of the order of the spin correlation length Iq ~ C,- This qualitative estimation leads to the 
choice fi = TT^/^^, kq = vr/^. Such a choice of Q gives a correct description of two limits. 
If ^ tends to infinity (the spin system has long range order, T = 0), the coupling radius of 
polaron also tends to infinity, kq tends to zero and Qk ( |18|) is transformed to Qr from 



(^W{). It is clear that in the opposite limit of high temperature, when the spin correlation 
between neighboring sites is small, the system must be described by local spin polarons. 



As it is seen from (|2^ ) if ^ (~ Iq) tends to unity the value kq is close to tt. In this case 
a(R — R') = (5rr' and the operators Arj,? = 4,5,6 represent local spin polarons close 
to polarons A-£i^i,i = 1,2,3 (in particular, for kq = tt, A^q is the linear combination of 
^R,i; ^R,2, and Ar,3). As a result, we have a continuous description of the model from a 
finite temperature to zero one. 

The spin correlation length was determined taking the spin wave spectrum which 
was calculated for the frustrated Heisenberg model in the framework of the spherically 
symmetric approach ^^. If the spin system is not far from the Neel phase, ^ may be 



determined using the expansion of the spectrum close to the AFM point Q [^, which 



will be denoted by ,^q. Then we use the criterion kq = vt/^q. The explicit expressions of 
spectrum and ^q are given in Appendix A. The choice of the parameter kq is different 
when ^Q becomes small. The value of ^q depends strongly on the frustration parameter 
p. For p > 0.15 we have ^q < 2 -^ 3 (kq = ^/^q ^ 0.357r) and for such a frustration 
the real value of ^ can essentially differ from ^q. Indeed, for the frustrated case the 
spin-spin correlation functions Cr =< Sr(,Sr(,+r > have the following dependence on 

■TV. liTSx ' '^ySy 

Cn = mi(R)(-l)"^+"- + m2(R)[(-l)"^ + (-1)"^] (23) 

where mi(R) ^ m2(R) for p ^ 1 (Neel type phase) and mi(R) ^ m2(R) for p close to 
unity (stripe phase with gapless spectrum at points Qi = (±7r, 0), (0,±7r)). The values 
mi(R) and m2(R) are dictated by the spin spectrum gaps in the points Q and Qi. For 
intermediate values of frustration 0.15 < p < 0.55 these gaps of the spin spectrum are 
comparable, and in this p-interval it is meaningless to determine ^ by the expansion of 
the spectrum near the points Q or Qi. On the other hand, a local polaron structure 
is taken into account by the operators Ari, Ar2, ^r,3 and the operators Ar^4, Ars, 
AnQ are introduced as candidates to describe polaron states with large or intermediate 
radius Iq = 2-^3. This means that we should fix kq ~ 0.357r for p > 0.15 {kq = tt/Iq, Iq = 
2-^3) and take kq = tt/^q for p < 0.15. The calculated values of ^q and the adopted below 
relation between the values of the parameter kq and different frustrations are represented 
in Table 1. 

To find the spectrum of excitations in the framework of the operators ([TT|), ([171) we 
calculate the matrix elements for matrices D and K. The calculation of elements of the 
arrays D and K is usually lengthy and it involves the calculation of complex commutators 
for the operators A^^i and i?k,j- Some of these commutators cannot be expressed explicitly 
by two-site Green's functions and some approximations are needed. The matrix elements 
are considerably simplified for a one-hole approach which can be adopted due to treating 
the problem in the limit of low doping (we will consider doping n < 0.2). Then they 
are expressed in terms of two- and multi-site correlation functions of spin operators (in 
the cases considered below — two-, three-, four- and five-site correlation functions). By 
taking into account the spherical symmetry the three-site correlators can be reduced to 
two-site correlators, and five-site correlators can be reduced to four-site correlators. The 
four-site correlators can be reduced to the form 

^RlR2R3R4 = ((SriSrJ (SR3SR4)) , 



where all sites are supposed to be different. To calculate such correlators an approximation 



47 is used: 



VR1R2R3R4 — Cr^^Cr^^ + -Cr^.^Cr^^ + -CiJ^C'iJas- 

As a result the matrix elements are expressed over the static spin-spin structure factor Cq 
(Fourier transform of Cr). The explicit form of K and D matrix elements are reproduced 
in Appendix B. Typical sums in the matrix elements have the form: 



Ug 



qef2 qi,q2efi 



The expressions analogous to Ug and V^i contain one and two summations on q G il. 
Each sum on q is proportional to Q/n"^ which is a small parameter in our approximation: 
as mentioned above Q/ir'^ = Kq/^^ — O-l (^^^ Table 1). This smallness allows to justify 
the approximation which consists in neglecting some terms proportional to Q^. 



IV. RESULTS AND DISCUSSION 

After solving the system ([T0|) for the chosen set of operators Ar j the resulting Green's 
functions have the form 

6 ^W (lA 

G,,{u;,k)=Y: ^"^^ ' , ^,J = l-6. (25) 



:{uj-ei{k)' 



According to ([T^) the residues zJ-^-.(k)^ zJ2)(^) determine the spectral weights n}^^^(k) of 
bare oxygen holes with fixed spin a and momentum k in the quasiparticle state |k, a, /) 
of the quasiparticle band ^^(k): 

Let us remind that the bare hole spectral weight satisfies the sum rule J2i iT-hai^) ~ ^ and 
the maximum number of holes per cell is equal to four despite the presence of six bands. 
This means that in this model the Luttinger theorem is not fulfilled. 

We present the results of the spin-polaron spectrum calculation for a temperature 
T = 0.2J. Let us mention that the spectrum has a weak temperature dependence up to a 
temperature T ~ 0.4J. Since we are interested in relatively small doping values x < 0.2 
we represent results for the two lowest bands / = 1, 2 only. In the high doping regime our 
approach may be insufficient since it is based on the one hole approximation and does not 
take into account the interaction between polarons. 

In Fig. 1 the spectrum £:i(k) and the bare hole spectral weight n)jg.(k) of the lowest 
band for different values of the frustration parameter, p = 0.05, 0.1, 0.13, 0.15, 0.2 and 

9 



0.25, are presented by contour plots taking r as the unit of energy. As mentioned already, 
we suppose some qualitative equivalence between doping x and frustration p. 



Fig. 1 
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Fig.l: Spectrum £:i(k) (figures a) and bare hole spectral weight n|j^(k) (figures b) for 
the lowest band presented as contour plots (in units r = 1) for different values of the 
frustration parameter p (0.05, 0.10, 0.13, 0.15, 0.20 and 0.25) corresponding to Figs. 
1.1 ^ 1.6. The symmetry points in the first quadrant of the BZ are denoted as F = 
(0,0), X = {(7r,0),(0,7r)}, M = (7r,7r). 
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Fig. 2: Electronic spectrum of tlie two lowest hole bands for the same frustration param- 
eters as in Fig. 1 along high symmetry lines of the BZ. Zero energy corresponds to the 
Fermi level and the value of r has been set to 0.4 eV. The parts of both bands which 
have considerable spectral weight n)^^ (k) > 0.05 are plotted by solid lines, the parts with 
small weight by dotted lines. 



The relation adopted below a; <^ p is given in Table 1. Of course this relation is a 
phenomenological one and it may need some rescaling. Nevertheless, the main conclusions 
are preserved if we admit the realistic assumption that the small and optimal dopings 
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correspond to p ~ -=- 0.1 and p ~ 0.1 -=- 0.25, accordingly. The contours decorated by full 
circles correspond to the Fermi surfaces calculated for the hole dopings x from Table 1. 
The diameter of the circles is proportional to the spectral weight n\i^l(k) along the FS. 
In Fig. 2 we show the electronic spectrum for the two lowest bands along high symmetry 
lines for the value of r = 0.4 eV as adopted in Sec. 2. The zero of energy corresponds 
to the Fermi level. These parts of both bands which have considerable spectral weight 
n\'^ (k) are plotted by solid lines, those parts with small n),^ (k) < 0.05 by dotted lines. 

Let us first discuss the dielectric (insulator) or heavily underdoped regimes which are 
represented by Figs. 1.1, 1.2 and 2.1, 2.2. In Figs. 1.1. a and 1.2. a the minimum of £:i(k) 
is close to the A^-point and the spectrum is rather isotropic near the band bottom. The 
dispersions along the directions F — M, and T — X — M reproduce the ARPES results 
(compare for example Fig. 2.1 with the dispersion in Fig. 3 of Ref. 0). The band width of 
the first band Wi ~ (4.4 — 3.3)r = 0.44 eV reflects also the experimental ARPES results 
Wi ^ 0.2 eV for BiaSrsCaCuOe+o.s 0, VTi ~ 0.3 eV for CaaCuOaCla ||, or VTi ^ 0.35 eV 
for La2Cu04 ^. Some of the uncertainty in these experimental values might be related 
to the drop of spectral weight near the F-point. 

The most important result in the heavily underdoped regime is the sudden drop of 
spectral weight in the lowest band as k moves from A^ to M, see Figs. 1.1. b and 1.2.b. The 
position of the k-line where this drop occurs is close to that one which gives the remnant 
FS in the ARPES experiment 0. The spin-polaron spectrum in Figs. 1.1. a and 1.2. a (2.1, 
2.2) has a symmetry close to the symmetry of the magnetic BZ, but the spectral weight 
n|j^(k) has the symmetry of the initial BZ, see Figs. 1.1. b and 1.2.b. 

For the case p = 0.25, x = 0.19, Figs. 1.6., 2.6, which we treat as close to optimal 
doping, a large FS centered at (vr, vr) is seen. It has a form consistent to the Luttinger 
theorem but with full filling (1 + x), not x. The spectral weight map from Fig.l.6.b 
demonstrates that such a large FS takes place due to small bare hole spectral weight at the 
spin-polaron k-states under the FS. The average value oin\^l(k) is close to 0.17 -C 1. Let 
us note that the average value of n^^(k) for k below the FS would be approximately two 
times greater if we should take the local polaron approach, i.e. restrict the basis operator 
set to (pTI). One would loose even a qualitative correspondence to the experimental FS. 

Let us compare the two lowest bands in the dielectric and the optimally doped cases 
along the A^ to F cut. In the former case, see Fig. 2.1, only the lowest hole polaron band 
is important (e.g. the upper electron one in Fig. 2.1), as the upper hole band has a small 
spectral weight. For the optimally doped case. Fig. 2.6, the second band is important 
as well. The correspondence to the ARPES experiment (see Fig. 3 in [Q]) is good if we 
assume that parts of the spectra are missed where the spectral weight is small (dashed 
lines in Fig. 2. 6) and if we treat the second band in the region close to F as an extension 
of the first one. Such a treatment is correct if we take into account the broadening of 
the bands. In our approach such a broadening would be equivalent to an additional 
splitting of bands due to additional operators (relative to the set (0), (|TB])). Such an 



improvement was investigated for the unfrustrated regular Kondo lattice model |^. In 
the optimally doped case {p = 0.25) the effective bandwidth W is W = e2,max— ^f ~ 0.55 
eV, see Fig. 2.6, {W ^ 0.38 eV for nearly optimally doped BiaSraCaCuOe+i |,||,||]). 
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The comparison with the dielectric state demonstrates some band narrowing with doping 
reduction as it is also given by ARPES 0. 

An important feature of the spectrum is the presence of a second band £2(k) well along 
the A^ to M cut, see Fig. 2.6. This branch resembles the main branch along the A^ to F 
cut if we shift it in the k-space by Q. Relative to the main branch along the A^ to M cut 
the band e2(k) is also shifted in energy by 0.1 -=- 0.2 eV. This branch is present as well in 
the underdoped case, see Fig. 2.4. An analogous branch (but without energy shift) was 
recently reported for underdoped Bi2212 |ll|] and it was treated as the shadow band first 
observed in Ref. [|I^. For the optimally doped compounds the presence of this band can 
be important for the plasma frequency at about ujp = 1 eV if one takes into account the 
corresponding interband transitions. Let us mention that in the ordinary band theory 
such transitions start at high energy Ri 1.2 eV and lead only to a constant contribution 



in the static dielectric constant 50 



Our calculations reproduce also the extended saddle point which is stretched in the 
direction from X to F close to the FS (see Fig. 1.6. a). Such an extended saddle point is 
observed in the photoemission data of optimally doped compounds and there has been a 
lot of speculations that it is important to the physics of the superconductors 0. Note 
that for this flat band region (close to (27r/3, 0)-point) the bare carrier spectral weight 
rif^lek.) is not small and is close to 0.22, see Fig.l.6.b . 

The underdoped case is represented in Figs. 1.3, 2.3 with p = 0.15, x = 0.06 and in 
Figs. 1.4, 2.4 {p = 0.15, x = 0.11). As we move from doping 0.02 to doping 0.11 (compare 
in consecutive order Figs. 1.2.b, 1.3.b and 1.4.b), the topology of lines with equal spectral 
weight changes dramatically. The remnant FS disappears, there appears finite spectral 
weight n/j^(k) close to the M point, the curvatures of lines change sign in a broad region 
and the spectral weight along the cut X — (vr, 7r/2) decreases. Between X and (vr, n/2) 
the spectral weight is transferred to the second band. Since that is also the region where 
the FS crosses the X-M-line with increasing doping, the nj^l(k) spectral weight along the 
FS is strongly changed. Such a strong decrease of spectral weight as k moves along the 
FS from the central region around N = (7r/2, tt/2) to points close to X = (tt, 0) can be 
treated as the opening of a high energy pseudogap near the point (vr, 0). 

This pseudogap 6 is often determined as 5 = s{7i, 0) — Sf (see 0). In our underdoped 
case {p = 0.15, x = 0.11) n|j^(7r, 0) is small and the first band should be only weakly 
present at the X point in ARPES. So, we must treat e2(7r, 0) as e{7i, 0) in the expression 
for the gap 6 (see Table 1). Then for the adopted above value r = 0.4 eV the value 
of the pseudogap is equal to 0.19 eV in accordance with ARPES results (0.1 -^ 0.2 eV 
P,p!2|-p!^). The dependence of the spectral weight n\^l{Tc,0) on doping given in Table 

1 is non monotonic. At first, a strong drop of n)j^(7r, 0) takes place in the underdoped 
region, but the spectral weight is sufficiently large in the optimally doped regime. Such a 
behavior reflects the ARPES scenario for a decreasing pseudogap as the doping increases 
from 0.1 to 0.2. 

The evolution of the lowest band dispersion shows that in the just mentioned transition 
region also the formation of the extended saddle point takes place close to the (27r/3, 0)- 
point. As the doping increases above x = 0.11 this fiat band region is preserved. 
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(1) 



For any doping, the particle spectral weight n^^(k) for k close to the A^ point is large 

and nearly constant nj^lCk) = 0.22 ~ 0.25, see Table 1. And indeed, in this k-region the 
ARPES experiments demonstrate the position of the FS surface with a well pronounced 
quasiparticle peak. Accordingly, the second band in this k-region gives a small value 
n\[^(k) < 0.03 (a small incoherent part) at any doping. 



Table 1: Frustration parameter p, doping x, spin correlation length ^q and cut-off 
momentum k for which calculations have been done. The derived energies (^^^-'(Tr, 0), 
e^'^^TT, 0), and ep) are given in units of r. n]^'^ (k) denotes the spectral weight of the first 
or second band, respectively. 



p 


Kq 


^Q 


X 


<i(vr,0) 


nS(vr,0) 


n«(7r/2,7r/2) 


£«(7r,0) 


£(2)(vr,0) 


ep 


0.05 


O.OSvr 


12 


0.00 


0.23 


0.02 


0.23 


-4.06 




(-4.42) 


0.10 


O.ITtt 


6 


0.02 


0.25 


0.00 


0.26 


-4.09 




-4.38 


0.13 


0.307r 


3.4 


0.06 


0.13 


0.11 


0.27 


-4.21 


-4.04 


-4.34 


0.15 


0.357r 


<3 


0.11 


0.06 


0.18 


0.23 


-4.38 


-4.03 


-4.40 


0.2 


0.357r 


<2 


0.14 


0.08 


0.16 


0.24 


-4.35 


-4.08 


-4.40 


0.25 


0.357r 


< 1 


0.19 


0.10 


0.13 


0.25 


-4.28 


-4.10 


-4.36 



V. CONCLUSION 



We presented a semi-phenomenological approach to describe the evolution of Fermi 
surface and electronic structure with doping. The approach is based on the spin-fermion 
model for the Cu02 plane and simulates the doping by a frustration term in the spin 
Hamiltonian. It shows a strong alteration of the electronic structure which is due to a 
simple and natural mechanism: doping leads to frustration in the spin subsystem and 
to the variation of spin-spin correlation functions, and correspondingly to the non-rigid 
behavior of the spectrum. 

Our theory explains the isotropic band bottom, the large energy difference between A^ 
and X [^0 and the remnant FS [^ of the undoped compounds. In the optimally doped 
case it shows a flat band region (an extended saddle point) between (vr/2, 0) and (vr, 0), a 
large FS |^|-0 and it shows signatures of a shadow FS ||10| , pl1| . We find in the underdoped 
region a rather strong change of dispersion and spectral weight at the FS near X. This 
finding can also be interpreted in terms of a pseudogap [[12|-[T^ . 

Our theoretical results would suggest for very low doping a FS in the form of small 
hole pockets (see e.g. Fig. 1.3. a). However, the spectral weights at those parts of the 
FS which are close to M are rather small. This might be a reason which hinders its 
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observation in ARPES. In underdoped samples we find a transition to a large FS but 
with very small spectral weight near X = (7r,0). Accordingly, we would interpret the 
Fermi surface arcs |jl5|] as those parts of the FS with large spectral weight. We find the 



parts of the FS with large spectral weight near to A^ to be remarkable stable with doping 
despite the tremendous changes in other parts of the BZ. 

There are several possible improvements of the approach presented here. A straight- 
forward way consists in introducing of several non overlapping domains J^i, Q2, ^3, ■■■^i 
instead of the single domain Q, (see also |^). As a result such an extension of the opera- 
tor set should lead to a better description of the incoherent part of the spectral function. 
It might also lead to a better description of the shadow FS and to a more abrupt change 
of the spectral density along the FS in the underdoped case, i.e. a better description of 
the pseudogap. 
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Appendix A 

Let us present the explicit form of Green's function (^^(k, uj) and spin- wave spectrum ci;(k) 
for three degenerate modes in the framework of the spherically symmetric approach ||45|| : 



G,ik,co)=<S^\S\>= /^^l,, , a = x,y,z (27) 



F(k) = -2 [Ji^,C,(l - 7,(k)) + J2ZdC,il - 7d(k))] (2J 



Uj\k) =|[ (1 - 79) {jlJ2Kgd + ll{Zg{Z, - l)Cga, + \z g + Kgg] 

+ (1 - 7d) [jiJ2Kgd + Ji{zd{zd - l)Cda^ + Izd + Kd, 
- (1 - 7g) -fdJiJ2ZgZdCga^ - (1 - 7d) -fgJiJ2ZgZdCda^ 



(29) 



317^-92 diyt-d2 

^gd = X! OirCr ; Kgg = ^ arCr ] Kdd = JI ^^^^ " 
r=g+d r=gi+g2 r=di+d2 

Here Zg, Zd are a number of first and second nearest neighbors on square lattice and a^ 
are vertex corrections. 
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ag = 


a^;ad = a^;ar = a^ if 


r > d; 






aj— 1 


= Ra = 0.863; 








«3 = 


(i-p)"2 +P"i; 








Jl = 


(1 -p)J,J2 = pJ; 








Kg, 


= 8a,Cg + 8a,Cf,Kgg = 


= 4^2 Csg 


+ 8a 


A 


Kdd 


= 8a,C2g + ^a^C2d, 









The values of a^, Cg, Cd-, C2g, Cgd, C2d are calculated self-consistently for each values of p 



and T [^ . The equations p7| . |28| . |29| ) give an expression for the spin-spin structure factor 
C(k): 



C(k) = A{1 + e"''/^)/(e"''/^ - 1), A = F(k)/2a;k. 



(30) 



The correlation length ^q (relative to the Neel type phase) is determined by power ex- 
pansion of the Green's function (^Tf ) on the value q=Q — katci; = ^^ : 



2 „2^ 



G(q,0) = G(Q,0)/(l + e4q 



(31) 



Appendix B 

Here we present the explicit form of K and D matrix elements. The following notations 
are used below. 

A,(k) = ({[Ak,„ (f + J + /i)], Al^}) = r% + Jij\f + J2j\f - hkr, 
Kij , fij , Jjj- , Jj\,- , /ijj are symmetric matrices. 

7g = 0.5(cos/cs + cos/cy); 7^ = cosk^ cos ky] 72c, = 0.5(cos2A;x + cos2A;y); 
The nonzero matrix elements are as follows: 

K — matrix 

Kii = 0.5(1 -|- cos kx); 

KiQ = 0.5m(1 -|- coskx); 

K22 = 0.5(1 + cosky); 

K2Q = 0.5^(1 -|- cos ky); 

Ks3 = 3/A + Cg^g; 

K34 = 0.5{u + Ug COS kx); 

K35 = 0.5(m + Ug cos ky); 

K36 = -u + jgCg{vg - v); 

K44 = 0.5(m -|- Ug cos kx); 

K4Q = —uv + 0.5w + 0.5 cos kx{—VgU — vUg + Wg); 

i^55 = 0.5(m -|- Ug cos ky); 

K5Q = —uv + 0.5w + 0.5 cos ky{—VgU — vUg + Wg); 

KqQ = (3/4)m + 2UV-W+ -fglugCg + m2 - (1/3)^2 - 4:Cg{{l/3)UgV + UVg) 

+2vvgCg + 2vlCl + Cgv\l + (2/3)C,)) - (2/3)7,W^M; 
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r— matrix 

Til = 0.5(1 + cosA;^)^; 

ri2 = 0.5 + 73 + 0.57^; 

ri3 = (3/2 + 2C373)(l + cosA;,); 

ri4 = (1 + cos kx){u + Ug cos kx)] 

ri5 = u{l + coskx) +Ug{cosky + 7d); 

Tie = (1 + cos A;j,)(m(-1 + 73) + 2Cg-fg{vg - v)); 

T22 = 0.5(1 + cos A; 



\2. 



y) 



T2Z = {l+COSky){'i/2 + 2Cg^g)- 

T24 = u{l + COS ky) + Ug{cos kx + 7d); 

T25 = (1 + cosky){u + UgCosky); 

T26 = (1 + cosA;j/)(m(-1 + 73) + 2C373(t;3 -t;)); 

Tss = -9/8 - ACg-fg + Cg + 0.5^2^725 + Cdid, 

Ts4 = — (3/4)-u — 0.5vCg — Ug cos k^ + (cos^ A;^ — 0.5) (f 236*3 + 0.5u2g) 

+ lg{Ug + '^VgCg - 2VgCg COS t^ " 2V C g) + 7^(0. 5^^ + V dC g) ] 

T35 = — (3/4)'u — 0.5f C3 — Ug COS /Cj^ + (cos^ ky — 0. 5) (f 23(73 + 0.5^23) 

+ lg{Ug + '^VgCg - 2VgCg COS ky " 2VC3) + 7^(0. 5^^ + V dC g) ] 

r36 = (9/2)m - (1/4)^3 + 73(3m + 4t;C3 + 3^3 - 4t;3C3) 

+l2g{CgU + 0.5t;23C23 - (2/3)t;23C3C23 - 0.5t;C23 
-{2/?,)vCgC2g - {l/^)UgC2g ' 2v gC] + U2gCg) 

+7d(2C3M + VdCd - vCd - iA/3)CgCdivd + v)- i2/3)ugCd - AvgC^ + 21.^^3); 
r44 = (3/4)-u — VgUg — 2uv + w + 0.5^^"^^ + cos A;^('U3 — 2f ^3 — 2f3-u + 21^3) 

+ (cos2 kx - 0.5)((1/2)m29 - 2t;3M3 + wf)- 
T45 = 0.5-u + 2^g{-vug - VgU + w g + .5u g) - 2uv + W + 7(^(0. 5^^ - 2-U3t;3 + wj^'^); 

TiQ = l.5UV + 1.5m - 0.75W + 0.5(mC3 + CgV"^ + CgVVg + CgV"^ - i8/3)UgVgCg 

+ {2/3)uUg - {8/3)uvCg + {8/3)vVgC^g) - (l/3)W^i^"^ 

+ -fg{-VgU - VUg + Wg + 2UgCg + 2V^ C g + {A/3)v^C^ + AVV gC g + AV^gC^ 

+2u' - {2/3)ul - (SmugvCg - SuVgCg - {A/3)W^^^) 

+ COS kx{—Wg + 1.5-U3 + UVg + UgV) 

+ {2cos'^kx — l)(0.25w2g — 0.25-ut;23 — 0.25-U-U29 + 0.5C3-U29 + 0.5C3t;t;23 
+ClvgV2g - {l/3)CgU2gV - C3Mt;23 - (1/6)m3M29 + O.5C3OT3 + 0.5C3t;3 

+ {l/3)ClvVg - {A/3)CgUgVg + Q .bUU g " (1/3)^^^")) 

+27rf(0.25wrf - Q.2huvd - 0.25vud + 0.5CgUd + 0.5^3^^ + C^VgVd - {l/3)CgUdV 

-CgUVd - {llQ)UgUd + Q.hCgVVg + 0.5C3t;| + {ll3)ClvVg - {A/3)CgUgVg 

+Q.buug-{l/3)wi% 

T55 = (3/4)-U — VgUg — 2UV + W + 0.5W^^) + cos ky{Ug — 2VU2g — 2VgU + 2Wg) 

+ (cos2 ky - 0.5)((1/2)m23 - 2t;3M3 + M;f ); 

ne = 1.5UV + 1.5m - 0.75W + 0.5{uCg + CgV^ + CgVVg + CgV^ - {8/3)UgVgCg 
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+ {2/?>)uUg - (S/?>)UVCg + {S/?>)VVgCl) - {l / 'i)W[^^ 

+ lg{-VgU - VUg +Wg + lUgCg + 2V^ C g + {^ j 'i)V^ C] + ^VV gC g + 4t;|C| 

+2«2 - (2/3)«2 _ (8/3)«,t;C, - %uVgCg - (4/3)W-W) 

+ cos ky{—Wg + 1.5-Ug + UVg + UgV) 

+ (2 COS^ ky — l)(0.25W2g — 0. 25^^23 ~ 0.25vU2g + 0.5Cg-U2g + 0.5CgVV2g 

+ CgVgV2g - {l/3)CgU2gV - CgUV2g - {l/6)UgU2g + 0.5CgVVg + 0.5CgVg 
+ il/3)C^gVVg - iA/3)CgUgVg + .5UU g " (l/S)^^^")) 

+27rf(0.25wrf - 0.25uvd - 0.25vud + O.bCgUd + O.bCgVVa + C^VgVd - {l/3)CgUdV 

-CgUVd - {l/Q)UgUd + Q.bCgVVg + ^ -bC gV] + {l/3)ClvVg " {A/3)CgUgVg 

+^.buug-{l/3)wi"^] 



Tee = -'^vu - (9/8)m + (9/2)u; - (l/4)u;« + ^(Cg + {2/3)ug - (8/3)C 



t^ 



WCg + t;2C, + t;t;,(C, + (8/3)C2) + t;,(0.5M, - (8/3)C,tz,) - (2/3) W^i^'^ 
+73(mg(-6 + (16/3)^3) + Vgu{-Q + 16CJ - ^vVgCg - v^{ACg + (8/3)^2) - Au'' 
+ (4/3)w2 - Sv^Cl - 4«,C, + 6t^, + (8/3)iyW) 

+725(0.5^25^2, + t't'25(C25 + {Al3)CgC2g) + t;2(0.5C25 + (4/3)C,C25 + (VS)^^!,) 
+^l,(^l, + (8/3)C,C2,) + t;(-(2/3)«2,C72, + (2/3)«,C2, - 2u2gCg - 2uCg) 
-V2g{2uC2g + {2/3)ugC2g + (2/3)m2c;C3 + ^uC g) + O.Sm^ - (1/6)M4 - {lA/3)vlCl 
+vVg{ACl - {A/3)CgC2g) + VgV2g{S/3){Cl - CgC2g) + Vg{{2/3)ugC2g + (2/3)«2,C, 

+2mC,) + 2w2gCg - {l/3)wf)C2g - (l/3)W^i;)) 

+27,(0.5M,Crf + t;t;,(C, + (4/3)C,C,) + v\<d.hCd + (4/3)C,C, + (1/3)C|) 
+^d(Cj + (8/3)C,C,) + v{-{2/3)udCd + (2/3)m,C, - 2iz,C, - 2uCg) 
-Vd{2uCd + (2/3)M,Crf + (2/3)m,C, + 2uCg) + O.Sm^ _ (i/g)^2 _ (i4/3)t;|C| 
+vVg{ACl - {A/3)CgCd) + t;,t;,(8/3)(C| - CgCd) + i;5((2/3)«,Crf + (2/3)m,C, 
+2uCg) + 2«;,C, - {l/3)wpCd - (l/3)W^i"^); 

j(i)_jnatrix 

J33 = -ACg + Cg73; 

J34 = 2Cg{yg - -u) + COS A;^Cg(0.5t; - 2vg + 0.5t;2g + Vd)] 

J35 = 2Cg(t;g - t;) + coskyCg{Q.bv - 2vg + 0.5t;2g + Vd); 

J36 = 2C3(t; - Vg) + 2^3 + ^g{Q.bCg{Vg -V)- 2Clv g - ^ .bU g - {2 j 3)U gC 2g 

-{Al3)ugCd + {Al3)vdCgCd + (2/3)t;2,C,C2, + {2l3)u2gCg + (4/3)MrfC,); 

J44 = 2Cg(t;g - t;) + Cg cos k^iS^.^v - 2t;g + ^.'bv2g + -^d); 

J46 = (4/3)C,K - ^) + (4/3)C^,w(t' - t',) + vliCg + (4/3)^2) + 2v''Cg + (4/3) W^.^/^ 
+WgCg(-3 - (4/3)CJ + 0.5cosA;^((8/3)m2 - (2/3)^^ - (2/3)^^23 " (4/3)MMd 
+ (8/3)Mt;,C, - [^l3)ugVgCg + (2/3)Mt;,C, + (2l3)u2gVgCg + (4/3)Mrft;,C, 
-(2/3)iz,t;C, - (2l3)ugV2gCg - {A/3)ugVdCg + Avpgi-0.5 - {2/3)Cg) 

+ AVVgCg - v'Cg " VV2gCg " 2VV dC g + VgV{0.5Cg + ( 2 / 3 ) C ^ ) 

+VgV2gCg{0.5 + {2/3)Cg) + 2vgVdCg{0.5 + ( 2 / 3 ) C, ) 
+ (2/3)C,«) + wW + 2wf - ACgWg) + (4/3)1^(3'^); 

J55 = 2Cg{Vg -V)+Cg COS ky{Q .^V - 2V g + Q .'bV2g + Vd)] 

J56 = (4/3)C,(w, - ti;) + {A/3)Cgu{v - Vg) + t;2(C, + (4/3)^^) + 2v^Cg + {A/3)W^i^ 
+vVgCg{-3 - {A/3)Cg) + . 5 COS A;j, ( ( 8 / 3 ) M ^ _ (2/3)^^ _ {2/3)uu2g - {Al3)uud 



+ {S/^)UVgCg - {S/?>)UgVgCg + {2 / ?>)UV gC g + [2 / ?>) U2gV gC g + {A/?>)UdVgCg 
-{2/2>)UgVCg - {2/?>)UgV2gCg " {A / 'i)U gV ^C g + AvlC g{-Q .h - (2/3)^3) 

+AvVgCg - v^Cg " vv2gC g " 2t;t;rfCg + Vgv{Q.^Cg + (2/3)^2) 
+^^,t'2,Cg(0.5 + (2/3)C,) + 2t;,t;rfC,(0.5 + (2/3)C,) 
+ (2/3)C,(t.f + t.« + 2wf) - ACgWg) + (4/3)1^(3'^); 
FO = -4MCg - (8/3)mm3 + SCgUv + (8/3)^^3 - {S/?,)wgCg + t;gMg((32/3)C3 - 4) + 2^^ 
+ {%/?,)vguCg + 3t;3C3 - 3t;C3 - %v^Cg - Qv^Cg - (8/3)t;|C2 + 2vVgCg - SC^vVg 

Fl = {2/3)u^ + UgCg - {10/3)uVgCg " 2M3t;C3 + 0.5ugV + 0.5Mt;g + v^Cg{0.75 + CJ 
+t;2Cg(0.25 + (5/3)6-3) + '^v^gCg - 0.5wg - {2/3)wg{C2g + 2Cd) + uvCg - uVgCg 

-{2/3)U2gVgCg + {2/3)UgVC2g - {2/3)UgVCg + {2/3)UVC2g + UgVgCg + {2 / 3)UV gC 2g 

-UgV^gCg - {2/3)uv2gC2g + {2/3yCg{Cg - C^g) + t;|C,( (5/3)C, - (7/3)C2,) 

-{S/3)VVgCl + VV23C<^((-5/3)C2g + CJ + AVgV2gCgC2g + 2(Mt;C<^ - UVgCg 

-{2/3)udVgCg + (2/3)Mgt;Crf - {2/3)ugvCg + (2/3)Mi;Crf + M3t;3Cg + {2/3)uVgCd 
-UgVdCg - {2/3)uvdCd + (2/3)t;2C,(C, - C,) + vlCg{{b/3)Cg - (7/3)C,) 

-(8/3)t;t;,C2 + i;^,C,((-5/3)C, + C,) + 4t;,t;,C,C,) + 2/3Cg{w% + 2^^^,^) 

-(2/3)W-W; 

F2 = SUgVgCg + 3UV gC g + {2/3)UVgC2g + {2/3)U2gVgCg " {2/3)UV2gC2g - UgV2gCg 

-3uvCg - {2/3)u2gvCg + {A/3)udVgCg - (4/3)Mt;,C, - 2ugVdCg - {A/3)UdVCg 
+ {A/3)uVgCd + ^.huvg - ^.huv + Cl{v2g + v + 2vd~ Avg) + w2((_ii/3)C| 
-0.25C, + (2/3)C2,C, + (4/3)C,C,) + vl{{-^Q/3)Cg^ 
-{l/3){CgC2g - {lA/3)CgCd - 0.25Cg) + t;t;,((52/3)C2 + 0.5Cg 

-{2/3)CgC2g - {A/3)CgCd) + VgV2g{{S/3)Cl + {A/3)CgC2g) 

+vv2g{{-h/3)Cl + CgC2g) + t;,t;.((8/3)C,C. + (16/3)^2) 
+vvd{{-lQ/3)Cl + 2C,C,) + (l/6)iyf ); 

j(2)_jjiatrix 

Jrfgs = -4Crf; 

Jd-ii = 2Cd{vd -v) + cos KCd{vg + Vf - 2vg)] 

Jd35 = 2Cd{vd -v) + cos kyCd{vg + Vf- 2vg); 

Jdse = 2Cd{v - Vd) + 2ud + {Al3)^g{Cd{CgVg + CfVf) - 2CgCdVg + Cg{ug + Uf) 

-ug{Cg + Cf)y, 

Jd44 = 2Cd{vd -v) + cos KiCdivg + Vf) - 2CdVg); 
Jd55 = 2Cd{vd -v) + cos ky{Cd{vg + Vf) - 2CdVg); 
Jd,4 = {A/3)uCd{v - Vd) + {A/3)Cd{wd - w) + 2v^Cd + vj{Cd + {A/3)Cj) + (4/3) W^i/^ 

-vVd{3Cd + (4/3)Cj) + cos A;,((4/3)C,(Mt;, - UgVd) + {A/3)ugUd - {A/3)CdW2g 

+VgVd{Cd - {A/3)Cj) + 2CdVg{v - Vd) - {2/3)u{ug + Uf) + {2/3)Cd{vd{ug + Uf) 

-Ud{vg + Vf)) + ((2/3)C| + 0.5Cd)vd{vg + Vf) - Cdv{vg + Vf) 

HV3)Ca{wfJ +wl^J) + {A/3)W^i^); 

Jd,, = {A/3)uCd{v - Vd) + {A/3)Cd{wd - w) + 2v^Cd + vj{Cd + (4/3)C|) + (4/3) W^j/^ 
-vVd{3Cd + {A/3)Cj) + cos ky{{A/3)Cd{uVg - UgVd) + {A/3)ugUd - {A/3)CdWg 
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+VgVd{Cd - (4/3)Cj) + 2CdVg{v - Vd) - (2/3)m(m3 + Uf) + (2/3)C,(t;rf(M, + Uf) 
-Udivg + v^)) + ((2/3)C| + ^.hCi)Vd{vg + t;/) - Cdv{vg + t;/) 

+ (2/3)C,(^fj + wf,) + (4/3)W^i3^)); 
Fo?0 = -^uCd - (8/3)««rf + 8Crf«t; + (8/3)u;Cd - (8/3)wrfCrf + t;dMd((32/3)Cd - 4) 

+ (8/3)t;rfMCrf + "iVdCd - "ivCd - Sv^C^ - GvjCd - {8/3)vjCj + Ivv^Cd - SCjvVd 

+2wll] + i8/3)iwt^ - W![P - W!^^); 
Fdl = -iA/3)wgiCg + Cf) + imCgiw^^] + wfj) + uv{ACd + (4/3)(C, + Cf)) 

+ {A/3)UgV{Cg + Cf) - {Sl3)UdVCg - AUVdCd + {A/3)UVg{Cg + Cf) 

-{A/3)u{vgCg + VfCf) - {A/3)vdCg{ug + u f) + (4/3)iz,t;,Crf + (8/3)M,t;,C, 
-(4/3)M,C,(t;, + t;/) - {2/3)ugCd{vg + t;/) + v\{^/3)CgCd - (4/3)C,(C, + C/)) 
-(32/3)^t;,C,C, + 2vCgCd{vg + t;j) - (10/3)t;C,(t;,C, + t;jC/) + t;,t;rf((20/3)C,C, 
-2C,(C, + C/)) + (8/3)C,t;,(r;,C, + VfCf) + (16/3)^,C,(t;,C, + t;;C;) 
-{S/3)vlCg{Cg + C/) + T;t;.((4/3)C,(C, + Cf) - (4/3)C,(C, + C;)); 
Fd2 = v\-{S/?>)CgCd + (4/3)^,(^3 + Cf)) + 16C,Crft;t;, - AC^uv + 2w,C,Crf 
+2t;t;jCdC/ - {lQ/3)vCdCg{vg + t;/) - (4/3)Cg(«, + Uf)v + (4/3)C,M,t;, 
+ (8/3)C,iZrft;, + (4/3)(C, + Cf)uvg - VgVd{2QCgCd + 2Cd{Cg + C/)) 
-(8/3)C,(C, + Cf)vl + ^;^;d((8/3)C,Crf - {A/3)Cg{Cg + C/)) + AC^uv^ 

+ VdCd{2Cg{Vg + t;/) - 2{VgCg + V/C/)) + (4/3)Cg(Mg + U f)V d - {2/3)CdUg{Vg + V f) 

-{A/3)u{CgVg + Cjt;j) - (4/3)C,M,(t;, + Vf) + (10/3)C,C,t;,(t;, + Vf) + (1/6)^^]-') 
+2CrfT;d(C3t;<^ + CfVf) + (8/3)C3W3(C3Wg + CfVf) - ACgCdVg + 2C3Crf(t;g + Vf); 
Jdee = FdO + -fg{Fdl + Fd2); 

h— matrix 

/ii2 = 2(0.5 + 7g + 0.57d); 
hi6 = 2m(0.5 + 7g + 0.57d; 

/i33 = 2(3/4 + 2C,79 + Cd7d); 
/i34 = 2(0. 5m + Mg7g + 0.5udjd); 

^35 = ^34; 

hsG = 2{-u + 2Cg7g(t;g - v) + Cdjd{vd - v)); 

/l45 = ^34; 

/i46 = 2(-Mt; + '^giwg - UgV - Vgu) + 0.5w + 0.57rf(wd - uvd - Udv)); 

h^e = /i46; 

hee = 2{{3/A)u + 2uv-w + ig{2CgUg + 2Cg{v^ + 2vgv) + {A/3)C^g{v^ + 3v^g) 

+2u^ - {2/3)ul - 8Cg{{l/3)ugV + uvg)) - (4/3)73W^i^) + ld{CdUd + Cd{v^ + 2vvd) 
+ (2/3)Cj(t;2 + 3^;^) - ACd{{l/3)udV + uvd) + u' - {l/3)ul) - (2/3)7.W^f ); 

Above the following notations are taken 



V = 

VI = 


- J_y-g«l. 


1 

N 

1 


E 

K 


g, 


d,2g,f: 


u = 


]v ^^«;' 
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ui = j^Ee'^'C,; l = g,d,2g,f; 

K 
Kl,K2 

«^/ = ^^E e— ia,_«,; l = g,d,2g; 

w\t = ^ T e-'-^'^e^^-'-C,,.^,- li,2= g,d,2g,f; 
m = 1 : Ua > 0; m = 2 : l^s < 0; m = 3 : I1I2 = 0; 

(m) (m) / T T \ 

Wi^^^ = ^ E e-*(''-''2)^e--^'C,C,_i; l = g,d,2g; 

-"^ Kl,K2iP 

W,['^ = ^ E^ e^^'a...^^^; 1 = g, d; 

W^z^-"^ = ^ ^E^ e^('^-''^)'a,_«2a2; 1 = g, d; 

W[s'^ = ^ E e— ^s(7/M-7/K-'^2))a,-K2C.2; l = g,d; 

Kl,K2 

Vr//^ = 1^ E (7K'^2) - tK/^i - t^2))C^,^n,C^,; 1 = g, d; 

W^/'^ = iT E E(Ee-^-^(C,+,C,_g-C,C,+,_g) 

Kl,K2 1 P 

_ ^ e-^('^i-'^2)Pg-«ig(7^C'^^j^^ + ^ e~^('^i-''2)Pg-«i(i-g)c'^C'^^j_^). 1 = g, d. 
p p 
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